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MICROWAVE EMISSION FROM SNOW AND GLACIER ICE 


I. INTRODUCTION ' 

Recent observations using the Electrically Scanned Microwave Radiometer 
(ESMR) on board the NIMBUS-5 satellite over the continental ice sheets of 
Greenland and Antarctica (Gloersen et al. , 1974) have been puzzling in that there 
is a complete lack of correlation between microwave brightness temperature and 
physical surface temperature in Greenland, whereas in Antarctica there is a 
loose correlation over wade areas. This has provided additional incentive for 
developing an applicable analytical model for explaining the microwave emission 
from snow and glacier ice. In addition, other experimental observations made 
on board aircraft sensor platforms over smaller glaciers, snow fields, and sea 
ice have been difficult to interpret in terms of existing models. 

A useful insight into the microwave emission from snow fields has been provided 
by a macroscopic volume scattering model (England, 1974), which was also 
developed for other solid dielectrics such as moist soil. England's treatment 
involves specifically a parameter called the volume scattering albedo, o) 0 , 
which is the ratio of the volume scattering coefficient to the total extinction 
coefficient. The extinction coefficient includes both the resistive losses and 
the scattering. The analysis involves specifying a value for gj 0 and then com- 
puting the brightness temperature or emissivity based on this single parameter. 
The model may be used to calculate the emissivity or the brightness tempera- 
ture of finite slabs of snow and ice with varying composition. 

In this paper, a microscopic approach will be taken in which it is assumed that 
the snow field or snow cover consists of closely packed scattering spheres 
which do not interact coherently. These assumptions then permit the utilization 
of the Mie scattering theory (Stratton, 1941) in computing w 0 and its several 
components. The treatment is similar in many respects to that used in calcu- 
lating microwave and optical properties of clouds (Gaut and Reifenstein 1971, 
Curran et al. , 1975). One important similarity is that the particle size is com- 
parable to the wavelength of the microwave radiation within the particle in both 
cases, since the particle sizes encountered in snow fields are sufficiently larger 
than cloud droplets to compensate for the lower index of refraction of the ice par- 
ticles. An important difference in detail of this treatment is that the imaginary 
part of the dielectric constant of ice is orders of magnitude smaller than for 
water droplets, grossly reducing the absoiption cross-section of the particle. 

Since the Mie scattering model provides many more reflecting surfaces than any 
reasonable multilayer (resulting from seaf onal snow accumulations) model or 
variable dielectric model would provide, ii; dominates the contribution to the 
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volume extinction coefficient at the microwave wavelengths and for subsurface 
snow particle sizes encountered. It will be shown that the penetration depths 
vary from 10 to 100 wavelengths depending on the snow and ice conditions 
encountered. Thus, the effects of snow cover dominate in most observations 
of snow-covered glaciers at the wavelength of 1, 55 cm used in the NIMBUS-5 
ESMR. 


n. ANALYTICAL APPROACH 


It shall be assumed that the snow and ice fields of interest consist mostly of 
spherical particles of water ice, with various values of the loss tangent at a 
given microwave wavelength. The exact theoretical study of the electromagnetic 
wave scattering by spherical particles was done by Mie (1908). By utilizing the 
Mie theory, the extinction and scattering cross-sections of spherical ice par- 
ticles can be calculated for all ranges of particle size and radiation wavelengths 
as a function of both the particle radius r and the complex index of refraction n. 

The extinction cross-section u EXT , scattering cross-section a SCA and absorption 
cross-section a ABS for particles of radius r and the complex index of refraction 
n can be written as: 


a EXT n- (2m + 1) ( l)Re[a m +b m ] 


2n 


( 1 ) 


m = 1 


°sca = -h (2m + 1) [|a m | 2 + |b m | 2 ] 


2 IT 


(2) 


m = 0 


and 


°ABS °EXT a SCA 


(3) 


where a m and b m are the complex Mie scattering coefficients which contain 
Bessel and Neumann functions of 1/2 integer order. That is, 
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where X is the free space wavelengtli, r is the particle radius, n is the complex 
index of refraction, J m denotes a Bessel function of order m and J m denotes the 
corresponding Neumann function. The extinction cross-section for spherical ice 
particles with refractive index n = 1, 78 + i 0. 0024, which corresponds to meas- 
ured values in the cm wavelength range of pure water ice near its melting 
point {Cumming, 1952), is calculated for several wavelengths and particle sizes. 
The results of the calculations in terms of the extinction coefficient (see Equa- 
tions (8) and (15) ), are shown in Figure 1. By using this relatively high value 
of hn(n), the importance of the scattering contribution to the extinction cross- 
section is emphasized, since it dominates even at this value of Im(n). The 
computation serves to illustrate the dependence of volume scattering on particle 
sizes in this wavelength range. 

Although the added complexity of calculating for various temperatures, snow 
density, and crystal size profiles is well within the capability of the computa- 
tional technique, a simplified model is constructed for the purpose of illustra- 
tion. A layer of uniformly sized spheres (the snow cover) of complex index n 
is situated on top of a perfectly absorbing medium (the glacier ice). The physi- 
cal temperature in the snow layer is assumed to be uniform and at 250° K for 
simplicity. The actual temperature of the underlying glacier ice, assumed to 
be 270° K, will be important only for very thin snow covers. Since the change 
in the dielectric constant from the ice to air would give less than 10% reflec- 
tivity and since the presence of the snow layer above both mitigates this abrupt 
dielectric change and introduces opacity between the ice and radiatively cold 
outer space, this model is reasonable except in the limit of va nis hing snow 


3 



WAVELENGTH (cm) 


Figure 1. Extinction coefficient as a function of microwave wavelength 

for several ice crystal radii 





thickness. The model is illustrated in Figure 2. Since the snow layer is 
assumed to be a temperature of 250°K, the appropriate index of refraction is 
n = 1. 78 + i0.0005 (Cumming, 1952). An additional simplification is that only 
the radiation emanating normal to the surface will be calculated, fix this case, 
the distinction between the two polarizations breaks down due to symmetry. 
Also, the ice particles are assumed to be randomly distributed and the inter- 
action between ice particles is not dependent on the path length. 



Figure 2. Sketch of the single layer scattering model 


The Rayleigh- Jeans approximation for the intensity of thermal radiation from a 
blackbody is applicable at microwave frequencies and at temperatures typical 
of the earth and its atmosphere; thus the radiative transfer equation (Chandra- 
sekhar, 1950) with axial symmetry, may be written as: 

2n 

T b (dj F(6,6 S ) Sin 0 s dd s + 7abs T(X) (7) 
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where T B ( 8 ) is radiance in the direction ( 8 ) expressed as equivalent blackbody 
temperature, T(X) is the thermodynamic temperature of the absorbing medium, 
and Text* T$ca* and Y^s are the extinction, the scattering, and the absorption 
coefficients respectively, defined as 


T i = 


No. 


( 8 ) 


where N is the number of particles per unit volume, F(0, 6 S ) is the scattering 
phase fmietion which is defined as the fraction of incident energy scattered into 
unit solid angle about a direction which makes an angle (H) with the 6 and 0 S : 


F(M S ) 


^-{|S,(®)I 2 + | S 2 «H)> t 2 } 
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where S (§) ) and S 2 ((H)) are scattering amplitude functions 
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a m and b m are fee same as defined in equations (4) and (5); n m and r m are func- 
tions of Legendre polynomials 


IT = P m (Cos 0 ) 

m m 7 
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Equation (7) is solved numerically by first noting that it may be expressed as: 


dT B (0) 

dX 





2tt 

T c (e s )F(0,0 s )Si n 0 s d0 s -7 SCA T B (0) 


(14) 


For a non- s catte ring medium, the left-hand side of the equation would be equated 
to zero; the right-hand side represents the angular distribution of radiation 
resulting from the scattering. To solve this equation, the snow layer is first 
divided into a number of smaller layers. An initial estimate is then made of 
the brightness temperature as a function of angle in each of these layers, insur- 
ing that die top and bottom layers obey the appropriate boundary conditions, 
namely, that the downwelling radiation is zero in the top layer and that the 
upwelling brightness temperature is 270° K in the bottom layer. The down- 
welling brightness temperature can then be calculated from the top layer down- 
ward using the estimated brightness for the scattering terms. The upwelling 
components are then similarly calculated from the bottom up. The process is 
repeated maintaining the boundary conditions until a satisfactory convergence 
is obtained, for which five passes was sufficient for all cases in this study. 

The method of solution is a variation of the Gauss-Seidal iteration (Hildebrande, 
1946). 


III. COMPUTATIONAL RESULTS 

Using the formalism developed in the preceding section, a set of brightness 
temperatures are calculated for wavelengths of 2. 81, 1, 55, and 0. 81 cm with 
ice particle layer thicknesses from 10 cm to 20 m and with particle radii 
ranging from 0.1 mm to 5 mm. The variation in the particle density, N(r) for 
different particle sizes is described as follows: 


N(r) = 



(15) 


This corresponds to mass densities typical of moderately packed snow (450 Kgm' 3 ). 
The calculated brightness temperatures are displayed in Figures 3, 4 and 5. 

In these figures, an asymptotic behavior i ? apparent when the thickness of the 
spherical particle layer becomes large compared with the wavelength, implying 
that little radiation emanates from depths greater than the penetration depth 
defined by such thickness. For the particle with radius 0.2 mm, typical of 
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BRIGHTNESS TEMPERATURE <°K) 








BRIGHTNESS TEMPERATURE (°K) 


BRIGHTNESS TEMPERATURE vs ICE CRYSTAL RADIUS 



Figure 4. Brightness temperature as a function of ice crystal radius 
for several values of the single layer depth and for a 
microwave wavelength of 1. 55 cm 
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BRIGHTNESS TEMPERATURE (°K) 


BRIGHTNESS TEMPERATURE vs ICE CRYSTAL RADIUS 



RADIUS (mm) 


Figure 5. Brightness temperature as a function of ice crystal radius 
for several values of the single layer depth and for a 
microwave wavelength of 0.81 cm 
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fresh snow, the asymptotic thickness is approximately 10 m for X = 2.81 cm, 

5 m for X = 1. 55 cm, and 2 m for X = 0. 8 cm. The asymptotic value of the 
brightness temperature of snow layers of increasing thickness depends only upon 
w 0 , the single scattering albedo in that layer. The asymptotic brightness 
temperatures calculated for 2. 81, 1, 55 and 0. 81 cm wavelength as a function of 
co o is plotted in Figure 6. 

The consequences of varying the index of refraction in the ice crystals has also 
been examined. The impurity level in the ice crystals is known to affect the 
real part of the index only slightly. On the other hand, the imaginary part of 
the index varies by many orders of magnitude with the introduction of impurities 
in the 260-273° K temperature range {Hoekstra and Cappillino, 1971). While it 
is approximately 0.0024 for pure ice near 0°C, it can be as high as 0.1 for 
first-year sea ice at the same temperature. Such variations are attributed to 
the lowering of the melting point of the ice and the existence of a two-phase 
system over a broader temperature range with the presence of metal salts. 

The variation of the imaginary part of the index of refraction leads to different 
values of cj 0 for a given wavelength X and particle radius r. This is illustrated 
in Figures 7, 8 and 9 for the three wavelength values 2. 81, 1, 55 and 0. 81 cm, 
respectively. These curves show that cu 0 is a strong function of the loss tangent 
which, in turn, depends on the physical temperature of the snow. This suggests 
that in order to determine the ice crystal size by remote sensing, multispectral 
measurements are needed. Curves such as those shown in Figure 10, in which 
co 0 has been eliminated as an intermediate parameter from the curves of Figures 
6-9, illustrates this requirement. Determination of the average particle radius 
from the curves of Figure 10 and from an a priori knowledge of the physical 
temperature range would be straightforward except for the complicating factors 
of particle size variation with depth (to be discussed later) and variation of 
penetration depth with microwave wavelength. 

Values of the azimuthaily integrated differential scattering cross-section, 

F(9, 6 S ) (see Equation (9) ), have been calculated and are shown in dimension- 
less form in Figure 11. The data are parameterized in terms of the particle 
size parameter a s 27rr/X and are divided by 2 ttX 2 so that these are valid for all 
wavelengths. Note that for a: - 0. 02, the forward and backward scattering are 
essentially equal and the scattering to the side is about half as large. As the 
particle size increases, the near forward scattering increases much more than 
the near backward scattering. (Note the curve for a = 1.12, Figure 11.) What is 
more important, however, is that the differential scattering cross-section 
increases rapidly as a increases, which accounts largely for the behavior of co Q 
as a function of the particle radius (Figures 7-9). 




11 



BRIGHTNESS TEMPERATURE 



EMISSl VITV 




CCT 


w vs. RADIUS FOR DIFFERENT LOSS 

° A = 2.81 CM 


0 . 002 ^ 


/ / / X^ .X* 0.015 

i // 

I / / / / ^ 

/ / / // / 

/ 1 / 

//A / / 

/ // 

' /// 

/7/x , 


- 0.050 /' 
X 
X 
X 

* — o.io 


.1 .2 .5 1 


RADIUS mm 


Figure 7 Volume scattering albedo as a function of ice cr y stel ra 
^ for several values of the imaginary part of the index of 

refraction and for a wavelength of 2. 81 cm 
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Figure 8. Volume scattering albedo as a function of ice crystal radius 
for several values of the imaginary part of the index of 
refraction and for a wavelength of 1. 55 cm 
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Figure 9. Volume scattering albedo as a function of ice crystal radius 
for several values of the imaginary part of the index of 
refraction and for a wavelength of 0. 81 cm 
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EMISSIVITY 


EMISSIVITY vs ICE CRYSTAL RADIUS 



RADIUS (mm) 


Figure 10. Emissivity as a function of ice crystal radius for three 
different values of the imaginary part of the index of 
refraction and for three different wavelengths 
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UNIT VOLUME DIFFERENTIAL SCATTERING CROSSSECTION 



270 ° 


Figure 11. Unit volume differential scattering cross section for three 
different values of the particle size parameter, or. The 
quantities shown in parentheses are the scale factors to be 
applied to the r-axis 


IV. OBSERVATIONS AND DISCUSSION 


An example of the spectral variation in the microwave emission of a snow- 
covered temperate glacier is shown in Tabl j 1. By neglecting the temperature 
gradient in the snow, the simplified model calculations can be used to determine 
the emissivity and consequently an average scattering particle diameter. The 
emissivity is taken as the ratio of the brigh ness temperature to the physical 
temperature of the surface, obtained both by remote measurements with an 
airborne 10 micrometer radiometer and frcm surface reports (Schmugge et al, 
1974). Assuming that the appropriate loss rangent for snow on a temperate 
glacier is 0.0014 (Im(n) = 0. 0024), the cun es of Figure 10 were used to obtain 
the particle radii listed also in Table 1. Note that the average radius increases 


17 



Table 1 


Microwave Brightness Temperatures of the 8.4 m Snow Cover on 
South Cascade Glacier (Station P-3, Elevation 2040 m, 1971 Data) 


X 

(cm) 

T 

a b 

(°K) 

e 

^physical = 267 K) 

^0 

r(mm) 

(for ny = 0.0005) 

0.81 

206 

0.77 

0.7 

0.15 

1.55 

194 

0.72 

0.8 

0.5 

2.81 

224 

0.84 

0.6 

0.5 

6.01 

238 

0.88 

0.5 

1.0 

11. 

256 

0.95 

0.3 

1-1.5 


with the wavelength. Since the optical depth is greater for the longer wave- 
lengths, the longer wavelength radiation is predominately from greater depths 
where the snow grain sizes are larger and the shorter wavelength radiation is 
from nearer the surface. These average radii are consistent with estimated 
grain sizes in the snow pack which range in radius from 0.15 mm near the top 
to 0. 5 mm at the bottom (Meier 1974). Although the model calculations pre- 
sented do not include the grain size, temperature, and density variations with 
depth, the simplified model does provide emissivities versus grain sizes which 
are consistent with the observations and at least qualitatively correct for wave- 
lengths up to 2.81 cm. For the longer wavelengths, part of the radiation 
emanates from the glacier itself, and curves such as those in Figures 3-5 apply. 

The computational results show that volume scattering from the individual 
grains is a dominant factor also in determining the emissivity of dry polar snow 
and fim, based on limited knowledge as to the particle sizes occurring in the 
snow cover. While this complicates the measurement of the near-surface 
temperature by microwave radiometry, it does provide a potential for measuring 
other parameters of glaciological interest, for example the extent of summer 
melting and snow accumulation rates. 

As pure water ice nears the melting point, the loss tangent increases rapidly, 
hi fact, for the thin liquid film that forms on the crystals, the loss tangent 
reaches a value the order of 0.4. Under these conditions, the scattering 
albedo, to 0 , is nearly zero (Figures 7-9) and the emissivity is nearly unity 
(Figure 6). This represents a plausible explanation of the sudden increase in 
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brightness temperature of snow observed near its melting point (Edgerton, 
et al, 1971; Gloersen et al 9 1974). This effect should enable the mapping of 
the annual variations of the dry snow zones versus percolation or soaked zones 
in Greenland, for example. 

A more detailed rationale is required to show the relation between snow accu- 
mulation rates and microwave brightness temperatures. It is known (Gow 1969, 
1971) that the crystal size in dry fim is a function of time, t, and temperature, 
T, 


D 2 = D 0 2 + k 0 (e~ E/RT )t 


(16) 


where D 2 is the cross-sectional area at time t, D 2 is the initial cross-sectional 
area, k 0 is a constant, E is the activation energy of the growth process per unit 
volume, and R is the gas constant. Also, the age, t, of a crystal at a particular 
depth, z, is a function of the snow load, o , and the mean accumulation rate, A, 
i.e. , 


q(z) 

A 


(17) 


o(Z), which is the integral over the density versus depth curve also depends 
weakly on T and A. To illustrate that the dependence of o(z) on A can be 
neglected, a comparison of ffi^.udheim (71.0° S, 10.8° W), where T = 256° K and 
A = 37 gm/cm 2 , with Wilkes (66.1°S, 110.6° W, where T = 254°K and A = 13 gm/ 
cm 2 , indicates that the snow load at Z = 25 m differs by less than 10% while A 
differs by a factor of 3 (Gow 1968, 1969). Therefore, for a given temperature, 
the crystal size depth profile depends on the accumulation rate, which is the 
parameter of interest: 


D 2 (z) 



(18) 


One must, of course, take into account the difference between crystals and 
grains, which are typically composed of one or two crystals (Gow 1969) and are 
a more likely measure of the scattering center size. However, grains and 
crystals have been observed to maintain a fairly constant mean size ratio of 
1.4:1 (Gow, 1969). 

hi general, the gradient of D 2 (Z) increases with increased T or decreased A, as 
can be seen from equation (18). The dependence on A and T is clearly indicated 
in Figure 12, where crystal size profiles obtained at three different field meas- 
urement sites are illustrated. 
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MEAN CRYSTAL RADIUS (MM) 



Figure 12. Ice crystal radius as a function of snow depth for several values 
of the accumulation rate, A and the average temperature, T. 
The values of the temperature are not shown with error limits; 
the minus refers to the smaller value of A and the + to the 
larger value. Note that the ice crystal radius variation with 
accumulation rate in each pair of curves bracketed runs counter 
to the trend that would be caused by temperature variation alone 
within that pair. 


Now for the range of grain sizes encountered, the corresponding portions of the 
emissivity curves of Figure 10 can be approximated by: 


e 


1 - e 0 D 2 


(19) 


where D 2 = nr 2 
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Since the curves of Figure 10 were calculated on the basis of uniform particle 
sizes, it is necessary to replace D 2 in Equation (19) by D 2 , which is an average 
of Equation (18) over Z, weighted by the radiative transfer properties. As A 
increases, both D 2 (Equation (18)) and D 2 will decrease and e (Equation 19) will 
increase. Therefore, emissivity and accumulation rates should be directly corre- 
lated over areas of constant temperature. 

We have in fact observed a remarkable correlation between the brightness tem- 
perature contours of the NIMBUS-5 ESMR (P. Gloersen et al., 1974) and the known 
accumulation rate patterns in Greenland (S. Mock, 1967) and Antarctica (C. Bull, 
1971), despite the fact that the spatial variations of the physical temperature are 
included. Since a 10% variation in emissivity is roughly equivalent to a 25°K vari- 
ation in temperature, it is not too surprising that emissivity variations dominate. 

In order to make a more direct comparison between emissivity and accumulation 
rate patterns, an approximate emissivity has been obtained, i.e. the ratio of 
brightness temperatures obtained at the wavelengths of 1.55 cm and 10 micro- 
meters. Such values cannot be obtained consistently because of persistent cloud 
cover, but were obtained on 11 January 1973 for Greenland (Figure 13) and Ant- 
arctica (Figure 14). Again, the areas with the highest emissivities are approxi- 
mately those with the highest accumulation rates. In making these comparisons, 
it is necessary to consider regions with similar physical temperatures due to its 
effect on grain growth rate. In particular, the emissivity in Greenland (Figure 13) 
decreases from the central summit (e = 0.85) to the Northern region (e= 0.70) 
which has lower accumulation. It is approximately constant (e = 0. 75) from 
Northwest of the summit toward Camp Century (see Figure 12) and increases 
from e - 0. 70 on the low accumulation area on the Southwest to e = 0. 85 on the 
very high accumulation area on the Southeast. In Antarctica (Figure 14), the 
emissivity is low (e = 0.70) over the low accumulation central East Antarctica 
plateau (vicinity of 80° S, 90° E). It is high (e = 0. 85) in the high accumulation 
area at the base of the Antarctic peninsula in contrast to a low value ( e = 0. 70 
to 0. 75) in the Ross Ice Shelf drainage basin part of West Antarctica (vicinity of 
82°S, 140°W). In both Southern Greenland and the Byrd station region of West 
Antarctica (80° S, 120° W), the T B contours are parallel to the ice divides which 
also mark regions of differing accumulation rates. 

In the case of the Lambert Glacier (72° S, 68 0 W) in Antarctica, the emissivity is 
also high, about 0. 85. It is known (W. Budd, 1974) that this glacier has extended 
areas of exposed blue ice due to vigorous wind action in the area; such areas 
would be expected on the basis of this model to have very low scattering albedo 
and hence high emissivity (see Figure 6). Thus, a high emissivity does not 
necessarily imply a high snow accumulation rate over a glacier. Therefore, 
some prior knowledge of the general characteristics of the snow field is required 
to interpret properly the observed brightness temperatures. 
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EMISSIVITY CONTOURS ON GREENLAND 
1 /11/73 


Figure 13. Contours of constant emissivity at a microwave wavelength of 1. 55 
cm on Greenland. The emissivity values were obtained by taking 
the ratio of the brightness temperatures obtained at wavelengths of 
1. 55 cm to those obtained at wavelengths of 10 micrometers from 
instruments on board the NIMBUS- 5 satellite. The data were 
obtained on a relatively cloud-free day (Jan 11, 1973). 
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EMISSIVITY CONTOURS ON THE ANTARCTIC CONTINENT 

1/11/73 


Figure 14. Contours of constant emissivity at a microwave wavelength of 

1.55 cm on Antarctica. The emissivity values were obtained by 
taking the ratio of the brightness temperatures obtained at wave- 
lengths of 1. 55 cm to those obtained at wavelengths of 10 micro- 
meters from instruments on board tie NIMBUS-5 satellite. The 
data were obtained on a relatively cloud-free day (Jan 11, 1973). 
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V. CONCLUSIONS 


A microscopic model, which considers the individual snow grains as the scatter- 
ing centers and employs Mie extinction and scattering coefficients, has been 
used to obtain a qualitative explanation of the low brightness temperatures 
observed over snow fields and snow-covered glaciers. More quantitative 
descriptions will be provided by this model when the vertical grain size distri- 
bution, physical temperature, and variation in loss tangent are taken into 
account. The present work has given essentially the same results as England 
(1974) in terms of the variation of brightness temperature with the scattering 
albedo to Q . The results differ in detail, however, because England assumed 
isotropic scattering in his model; Mie scattering, assumed here, is anisotropic. 
An additional result of the present analysis is the determination of the variation 
of oj 0 as a function of particle sizes, wavelength, and loss tangent. It has been 
shown that in the wavelength range from 0. 8 to 2. 8 cm, and for snow crystal 
sizes normally encountered, most of the microwave radiation emanates from a 
layer the order of 10 meters or less in thickness. An additional important 
result of this analysis is the explanation of the wet snow signature in terms of 
a wet surface on the snow crystals and the resulting sudden increase in the loss 
tangent of the seatterers which, in turn, gives rise to near-zero values for co 0 , 
the scattering albedo. Finally, it is concluded that it may be possible to deter- 
mine snow accumulation rates as well as near-surface temperatures by utilizing 
this model in conjunction with multispectral observations in the microwave 
region. 
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